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Description 

BACKGROUND OF THE INVENTION 

5 [0001] The present invention relates to an image processing method and an apparatus that prepare treatment plan- 
ning based on three-dimensional data obtained by making measurements on a patient using an X-ray computed tom- 
ography (CT) unit in treating diseases such as cancer with radiation. More specifically, the present invention is directed 
to an image processing method and apparatus adapted for radiotherapy treatment planning using digitally reconstruct- 
ed radiographs (DRRs), which are permeation images generated from three-dimensional data obtained from a \qvr\o- 

10 graphic image measurement apparatus. 

[0002] Treatment methods involving injection of beams of radiation such as x-rays and beams of protons onto a 
focus portion such as cancer are considered effective. To give a patient such a treatment, preliminary treatment planning 
must be prepared. A radiotherapy treatment is generally given under the following procedure. 
[0003] First, part of a patient body including an affected part is'measured using an apparatus such as an X-ray CT 

is unit. The affected part is specified from the measured data, and its position and size are grasped. Then, the isocenter 
is set to the affected part, and conditions such as the direction of irradiation, number of injections and range of irradiation 
are simulated and adjusted so that radiation can be focused onto the affected part as closely as possible. Then, based 
on the results of the simulation and adjustment, markings are made on the patient body. Thereafter, the patient is 
requested to go to the radiotherapy treatment unit, positioned on the unit in accordance with the simulated markings, 

20 and given the treatment. 

[0004] At the time of the aforementioned simulation, a method of preparing treatment planning using DRRs is avail- 
able. The DRR is a photographic image obtained by projecting onto a plane the pixels of data produced by a computed 
tomography unit such as an X-ray CT unit (these pixels will hereinafter be referred to as "voxels 1 *) using radially ex- 
panding rays that are irradiated from a radiation source. 

25 [0005] The treatment planning using DRRs provides the advantage that correct simulations can be made by calcu- 
lating the radiograph of an affected part based on the same paths as those of the radiation beam provided by the actual 
treatment unit. 

[0006] However, simple calculation of a radiograph imposes the problem that the transmitted body structure is hard 
to grasp. The following literature provide some solutions to this problem: (1) J. M. GALVIN, C. SIMS, "THE USE OF 

30 DIGITALLY RECONSTRUCTED RADIOGRAPHS FOR THREE-DIMENSIONAL TREATMENT PLANNING AND 
CT-SIMULATION," Int. J. Radiation Oncology Biol. Phys. Vol. 31, No. 4, pp. 935-942, (1995) and (2) JP-A-8-164217. 
[0007] Such literature disclose DRR generation methods that use a lookup table for making CT value conversions 
to highlight a bone area serving as a landmark when a radiograph is calculated. The CT value conversions are made 
using a bone window or the like. 

35 [0008] However, the conventional techniques are successful in improving the contrast of the bone area, but unsuc- 
cessful in generating a DRR having an area of interest highlighted more clearly. If clear highlighting is implemented, 
the user can prepare a more accurate treatment design. Further, while an exemplary method disclosed in JP-A- 
8-164217 involves an edge process to project an area such as a bone area clearly, the edge process is not specifically 
described. Thus, an effective edge process to enhance the contour using treatment data from an X-ray CT unit or the 

to like must be developed. 

[0009] Moreover, for treatment planning, measurements must be made with respect to the size of a target tumor and 
the size of an irradiation range. Since a DRR is a projected image formed by radially expanding rays, such measure- 
ments cannot be made on the DRR. 

[0010] Furthermore, it takes much time to generate a DRR. Therefore, for interactive data processing, a high-per- 
45 formance computer is required, and this has prevented interactive processing with a low-performance computer. 

[0011] EP 0 506 302 A1, US-4,737,921 and US-5,280,428 disclose to associate color information to three-dimen- 
sional image data on the basis of rays which are in parallel and are different from rays from a radiation source used 
to generate the three-dimensional image data. 

50 SUMMARY OF THE INVENTION 

[0012] To overcome the aforementioned problems, the present invention provides an image processing apparatus 
according to claim 1 for allocating color information so as to correspond to pixel values of three-dimensional data that 
is measured by an X-ray CT unit or the like, and allows a DRR to be generated by projecting the data using color 
55 information converted from the pixel values. In addition, the present invention provides an image processing method 
according to claim 9. 

[0013] In the present invention, a position for projecting a DRR may be on a flat plane which includes an isocenter 
and which is perpendicular to a line connecting the isocenter to a radiation source. In the present invention, the user 
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, may specify the number of pixels and resolution of a DRR, and the sampling interval of voxels on a ray connecting a 
pixel on the DRR and the radiation source. Further, intermediate images of low resolution may be displayed until a 
final DRR is displayed. 

[001 4] As described above, the image processing method according to the present invention allows the user to clearly 
5 display and easily grasp an area of interest such as a bone structure and the contour and shape of organs by means 
of an edge process when the user prepares treatment planning based on a DRR, which is a radiograph, consisting of 
three-dimensional data gathered from a computed tomography unit such as an X-ray CT unit before a radiotherapy 
treatment is given to a patient. 

[0015] Moreover, the image processing method according to the present invention allows the user to measure .data 
10 on a flat plane of a DRR including the isocenter of an area of utmost interest when data such as the size of a target 
tumor and the range of irradiation are measured. Furthermore, the method allows the user to select image quality and 
speed. Since the user can roughly grasp the condition of the area of interest before obtaining a final, high-resolution 
DRR, the user is allowed to process data interactively by, e.g., changing the direction of irradiation in order to generate 
the final DRR with improved resolution and increased number of' pixels. 

15 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0016] 

20 Fig. 1 is a flowchart of an image processing method, which is an embodiment of the present invention; 

Fig. 2 is a diagram showing the entire process of radiotherapy treatment planning; 

Fig. 3A is a diagram showing an exemplary unit for specifying a range of voxel values to which color information 
is allocated; 

Fig. 3B is a diagram showing an exemplary unit for allocating the color information; 
25 Figs. 4A to 4C are diagrams showing another exemplary unit for specifying a range of voxel values to which color 

information is allocated; 

Figs. 5A to 5D are diagrams showing an exemplary edge process based on a standard deviation using two-di- 
mensional data; 

Fig. 6 is a diagram showing an exemplary block for calculating a standard deviation for three-dimensional data; 
30 Fig. 7 is a flowchart of an image processing method for executing an edge process, which is an embodiment of 

the present invention; 

Fig. 8 is a flowchart of an image processing method for executing an edge process, which is another embodiment 
of the present invention; 

Fig. 9 is a diagram showing an exemplary projection made by setting a DRR plane at a position including an 
35 isocenter; 

Fig. 10 is a diagram showing an example in which the number of pixels and resolution of a DRR, and a sampling 

interval of voxels on a ray connecting a pixel on the DRR to a radiation source are specified; 

Figs. 11 A to 11 D are diagrams showing an example in which a high-quality DRR is displayed on a step-b^-step 

basis; 

40 Fig. 12 is a diagram showing an exemplary table for allocating color information; and 

Fig. 13 is a block diagram showing an image processing apparatus, which is an embodiment of the present inven- 
tion. 



45 
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DESCRIPTION OF THE PREFERRED EMBODIMENTS 



[0017] Embodiments of the present invention will now be described in detail with reference to the drawings. 
[0018] Fig. 1 is a flowchart showing a radiograph generation method, which is an embodiment of an image processing 
method to which the present invention is applied. The flowchart will be described with reference to the process of a 
radiotherapy treatment shown in Fig. 2. 
50 [001 9] The image of an area including an affected part is made using a tomographic image measurement apparatus 
1 such as an X-ray CT unit shown in Fig. 2 before a patient is subjected to a treatment, and thus three-dimensional 
image data is prepared (Step 101). Then, irradiation conditions necessary for treating the patient with a radiotherapy 
treatment unit 7 are set (step 102). The irradiation conditions include the position 3 of a radiation source, the isocenter 
position 4 onto which radiation is focused, the degree of expansion of a radiation beam 5 and so forth. The position of 
the radiation source is rotatable about the specified isocenter by rotating and moving a gantry 8 and a table 9 of the 
radiotherapy treatment unit 7. 

[0020] Successively, to highlight an area of interest, a color information conversion iabie corresponding to voxel data 
(hereinafter referred to as "voxel values") such as CT values is set (Step 103). The table contains color information 
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corresponding to a range of voxel values in terms of RGB ratio, as shown, e.g., in Fig. 12, where V fe the voxel value. 
The table shown in Fig. 12 can allocate values linearly within the range of voxel values. In this example, the voxel 
values ranging from -1000 to -100 are converted to the color information of gray having gradations ranging from 0 to 
1 28, the voxel values mainly indicating vascular areas and ranging from -99 to 99 are converted to the color information 
5 of red having a gradation of 200, and the voxel values mainly indicating bones and ranging from 100 to 1000 are 
converted to the color information of yellow having 0 to 255 gradations. 

[0021] An embodiment of a method by which a user specifies the color information table is shown in Figs. 3A and 
3B. First of all. as shown in Fig. 3A, the user specifies a range 11 for which color information is obtained. The range 
11 can be selected by, e.g., dragging both ends of a target range 11 using a mouse. Then, as shown in Fig. 3B, )he 

10 user prepares the color information table corresponding to the selected voxel value range. In this example, the user 
specification is such that colors G and B change as shown by a line 12 and that color R remains at 0 as shown by a 
line 13. These lines are specified by, e.g., dragging the ends of each line using the mouse. Through these operations, 
the color information table such as shown in Fig. 1 2 can be prepared automatically. Then, a DRR, which is a radiograph, 
is produced by sequentially calculating all the pixels on the DRR. 

is [0022] First, assuming a ray connecting a pixel on the DRR and the radiation source, voxel values on such ray are 
interpolated (Step 106 (since Steps 104 and 105 result in "NO" in the first operation, so that the process proceeds to 
step 106)). The interpolated voxel values are converted to color information by means of the color information table 
set in Step 103 (Step 107). Successively, pieces of the converted color information are added up to a memory for 
storing the pixel values of the DRR (Step 108). The pixel value storing memory is initially cleared, and the pieces of 

20 color information corresponding to the voxel values on the ray are sequentially added up. The addition is made on a 
color component (R, G, B) basis. This calculation is made for all the voxel values in the three-dimensional data present 
on the ray (Step 105). Upon completion of the calculation for all the pixels on the DRR (Step 104). the calculated DRR 
is displayed (Step 109). The DRR may, thereafter, be printed or stored in a hard disk, etc. Thus, this method, in which 
different areas are displayed in different colors, allow the user to discern bones, organs, tumors and the like with ease. 

25 [0023] According to the aforementioned method, different colors can be used to display different areas that are 
specified by respective voxel value ranges, and this allows the user to highlight areas of interest in preparing a treatment 
design. 

<Method of Specifying Voxel Value Range> 

30 

[0024] Another embodiment of a color information specifying method will be described with reference to Figs. 4A to 
4C. 

[0025] Using there-dimensional data obtained by a tomographic image measurement apparatus before treatment, 
a sliced image 21 is displayed. The position of the displayed cross-sectional image is indicated in the form of a line 23 
35 on a scanogram 22 so that the user can locate the position. The position of the sliced image 21 can be moved by 
dragging the line 23 with the mouse or the slicing position can be moved with a slider 24. 

[0026] Then, a single point in an area to be highlighted is specified on the sliced image 21 using a pointer 25. The 
voxel value V of the specified point is retrieved, and the retrieved value V is set as the central value of a voxel value 
range 26. so that the range 26 covering an interval of ± a from the central value V can thereafter be set for preparing 

40 a color information table. While a predetermined initial value is used as a, the interval a can be altered by dragging 
either the upper end or lower end of the color information setting range with the mouse. Once the voxel value range 
is specified in this way, the color information table is prepared based on the method described with reference to Figs. 
3A and 3B. This method allows the user to visually specify the voxel values of the area to be highlighted. 
[0027] Further, if the user specifies a voxel value range and displays the sliced image 21 using only the voxels 

45 contained in the specified voxel value range, the voxel value range specifying operation becomes easier. 

[0028] Still further, instead of specifying the area to be highlighted with a single point, such area may be specified 
by determining the voxel value range for color information specification based on the maximum and minimum values 
of the voxel values contained in a rectangular area. 

50 <Standard Deviation> 

[0029] A display method by which a structure such as a bone can be easily identified on a DR R by smoothly enhancing 
the contour of the structure such as a bone will be described with reference to Figs. 5A to 5D. To enhance the contour 
of an image, an edge process may be performed so that a voxel value at a position where there is an abrupt voxel 
55 value change is converted into a higher value than voxel values adjacent to such voxel value. An embodiment of the 
present invention, in which a DRR is prepared using edge-processed data, will be described below. Since techniques 
based on differential analyses produce data plenty of noise with highfrequency components emphasized, a technique 
based on a standard deviation is used in this embodiment to express desired voxel values. 
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, [0030] A standard deviation c can be given by equation 1 , assuming that the number of voxels within a block is n 
and that the density of a voxel i is Vj. 
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[Equation 1] 
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[0031] Figs. 5A to 5D show the results of an edge process that is performed with a block size of 3x3 using two- 
dimensional data. The 3x3 block is the minimum unit for performing the edge process, and the two-dimensional data 

20 is used for a concise description. Fig. 5C shows the results 32 of an edge process performed by calculating a 3x3-block- 
based standard deviation for an area surrounded by the thick line in two-dimensional data 31 shown in Fig. 5A. The 
graphs on the right show changes in the data. The graph 33 shows a change in the original data surrounded by the 
thick line and the graph 34 shows a change in the edge-processed data. The numerals given on the horizontal axes 
of the graphs 5B and 5D correspond to the numerals superscribed on the original data surrounded by the thick line. 

25 [0032] As seen from the graph 33, the edge in the original data surrounded by the thick line is assumed to be located 
between the second and third pixels from the left. It is verified that the contour is enhanced in the graph 34 that is edge- 
processed by means of a standard deviation with high values given between the second and third pixels as an edge. 
Further, there is another edge between the fifth and sixth pixels from the left in the original data. This edge, having a 
smaller density difference than the first edge, is expressed by smaller values in the edge-processed data 32. Another 

30 edge between the sixth and seventh pixels from the left in the original data, having a further smaller density difference, 
is expressed by further smaller values in the edge-processed graph. 

[0033] As in the aforementioned example, the edge-processed data is obtained by converting the original data as if 
the edge were present between voxels of the original data, i.e., between the second and third pixels of the original 
data. As a result, the contour is extended by one pixel. The conversion is performed so that the edge-processed voxels 
35 2 and 3 represent an edge. However, this technique provides smoothing effects, which reduces noise and hence allows 
the contour of the structure such as a bone to be enhanced satisfactorily. 

[0034] The edge process described in terms of two-dimensional data shown in Figs. 5A to 5D can be applied directly 
to the horizontal and vertical directions of three-dimensional data. 

[0035] Three-dimensional data is subjected to an edge process in the following way. When a standard deviation is 
40 calculated using a 3x3x3 block 41 around a voxel 42 that is subjected to an edge process, as shown in Fig. 6, satisfactory 
results can be obtained. It should be noted that, when the aforementioned block is used, the edge process cannot be 
performed on all the voxels on the outermost sides of the three-dimensional data. However, it is negligible since it is 
a small area in comparison with all data quantity. 

[0036] The edge-processed three-dimensional data is the edge-enhanced data in which portions of the original three- 
45 dimensional image with their densities drastically changed, i.e., edge portions are enhanced. 

[0037] Fig. 7 is a flowchart, which is an embodiment of the present invention. The flowchart presents a process using 
an edge process and a color information conversion table. The process steps identical to those shown in Fig. 1 are 
denoted by the same step numbers. 

[0038] In Fig. 7, a standard deviation a of the three-dimensional image data measured by a tomographic image 
50 measurement apparatus is calculated by the equation 1 , and a second version of three-dimensional data is obtained 
by expressing the original three-dimensional image data using a in Step 201 . In subsequent steps, the image process 
is performed using the second version of three-dimensional data. 

[0039] Step 201 may be performed immediately after Step 102 or immediately after Step 103. 
[0040] Fig. 8 shows another embodiment of a flowchart showing an edge process. The flowchart of Fig. 8 is distin- 
55 guished from that of Fig. 7 in that the color information conversion table is not used. Therefore, Steps 103 and 107 in 
the flowchart of Fig. 7 are eliminated, and Step 202 replaces Step 108. In Step 202, voxel values are added up to a 
memory for storing pixels on a projected image. Other steps are identical to those of Figs. 1 and 7. 
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<Position of DRR Plane> 

[0041] An embodiment of a method of projecting a DRR onto a position passing through the isocenter will be de- 
scribed with reference to Fig. 9. The DRR may be prepared from non-edge-processed voxel data or from edge-proc- 

5 essed voxel data. Fig. 9 is a cross section taken from a diagram showing the treatment planning of Fig. 2. The DRR 
is set so as to coincide with a plane 6 that includes the isocenter 4. The three-dimensional data 2 is projected onto the 
DRR plane 6 through the radiation beam 5 irom the radiation source 3. The three-dimensional data is added up and 
projected in both forward and backward directions indicated by arrows. That is, all the voxels in the three-dimensional 
data present on the ray connecting the radiation source to a pixel 51 on the DRR plane are added up and projected. 

10 The voxel values are projected after converted into color information using the color information conversion table shown 
in Fig. 12. However, when the edge-processed voxel data are used, the color information conversion table is not nec- 
essarily used. The DRR calculation method includes an image-order approach in which a DRR is calculated in the 
order of pixels or an object-order approach in which a DRR is calculated in the order of three dimensional data. 
[0042] According to the aforementioned embodiment, the projected position of the DRR is on a plane including the 

15 isocenter. Since voxels which are remoter from the radiation source than the projected area and considered invalid in 
ordinary projections are projected, all the three-dimensional data can be visualized as the DRR. If the user wishes to 
generate a DRR in which only an area of interest is focused, all the data may not necessarily be visualized by limiting 
the area for calculation. 

[0043] Further, the isocenter is usually set to the center of an affected part to be treated, and the range of irradiation 
20 is often set on a plane including the isocenter. Thus, in this embodiment, the irradiation range is displayed on a DRR 
by setting the pixel size of the DRR to coincide with the minimum value of the voxel size. As a result, the image of the 
affected part on the DRR becomes as big as the actual affected part. Hence, the DRR on the plane including the 
isocenter that is the center of an area of interest allows the user to make simulations easily. 

25 Specification of Resolutions 

[0044] An embodiment of a method of improving the quality and speed of image forming operation by specifying the 
number of pixels and resolution will be described with reference to Fig. 10. Similarly to Fig. 9, Fig. 10 shows how a 
DRR is generated in the form of a cross section. In this case also, the DRR can be generated by using a color information 

so conversion table and an edge process. 

[0045] First, the resolution (pixel size) 52 and number of pixels of a DRR are specified. A value equal to the resolution 
(voxel size) of the three-dimensional data is set as an initial value (the finest resolution value is set in case of that the 
voxel size has different lengths), and such number of pixels as to allow all the three-dimensional data as to be projected 
(such number can be obtained from the radiation source, isocenter position, number of three-dimensional data voxels, 

35 resolution, and the like) is specified. Then, a sampling interval 53 on the ray connecting the radiation source 3 to the 
pixel 51 is specified. The value equal to the resolution of the DRR is specified as an initial value. Using the aforemen- 
tioned initial values, the DRR corresponding to the quality of the three-dimensional data can be calculated. 
[0046] If the user wishes to save calculation time, the number of pixels of the DRR can be reduced, or the sampling 
interval between voxels on the ray connecting the radiation source to the pixel is increased. As a result, the DRR can 

40 be calculated faster. If the number of pixels of the DRR is reduced without changing the resolution, the projection range 
is reduced, but the image quality remains unchanged. On the other hand, if not only the number of pixels is reduced 
but also the resolution is reduced accordingly, the image quality is impaired, but the projection range remains un- 
changed. The above similarly applies to the sampling interval between voxels on the ray. Further, if the user wishes 
to improve the image quality, the resolution of the DRR may be increased, or the sampling interval between voxels on 

45 the ray may be reduced. 

[0047] As described above, the aforementioned embodiment allows calculation speed and image quality to be con- 
trolled in accordance with the usage of an image and the capacity of a computer. 

improvement of lnteractivity> 

50 

[0048] It takes time to generate a DRR. A solution to this problem is to improve interactivity by increasing the speed 
of generating a DRR in nominal terms. An embodiment of a method of implementing such high-speed operation will 
be described with reference Figs. 1 1 A to 1 1 D. Figs. 1 1 A to 1 1 D show a method in which calculations for generating a 
5x5-pixel-based DRR are made on four steps and the calculated image is displayed every time the calculation is made. 
55 in this case also, the DRR can be generated by using a color information conversion table and an edge process. 

[0049] Let us assume that a final DRR 61 is obtained through calculations. The numbers on the image are allocated 

io the respective pixels constituting the image, and indicate the order in which the pixels are calculated. 

[0050] First, pixels 1 to 9 are calculated every other pixel in both horizontal and vertical directions. A first-step DRR 



6 



EP 0 884 697 B1 



,62 is generated from the initially calculated nine pixels, and the generated DRR 62 is displayed" The first DRR is 
displayed with the skipped pixel values (e.g., 14,20, 10) represented as the calculated pixel values (e.g.. 1) (Fig. 11B). 
Then, pixels 10 to 13, each being located at a position one pixel moved in both horizontal and vertical directions from 
the initially calculated pixel, are calculated. A second-step DRR 63 is generated using the calculated thirteen pixels 

5 and displayed. In this case, the pixels on the left of the pixels 10 to 13 are displayed as the same values (Fig. 11C). 
Then, pixels 14 to 19, each being located at a position one pixel moved in the horizontal direction from the initially 
calculated pixel, are calculated to generate and display a third-step DRR 64 using the nineteen pixels (Fig. 1 1 D). Finally, 
the remaining pixels 20 to 25 are calculated to generate and display the final DRR 61 (Fig. 11 A). 
[0051] According to the aforementioned method, the user can grasp the outline of the DRR with the images 62, 63 

10 and 64 on the intermediate steps sequentially displayed before all the pixels of the final DRR 61 are completely cal- 
culated. Thus, interactivity can be improved. The total time required for completing all the calculations, although a small 
time to generate the intermediate images is included, is substantially the same as that required in the method in which 
a DRR is displayed after all the pixels are directly calculated. 

[0052] Fig. 13 is a block diagram showing an exemplary image processing apparatus that embodies the image 
15 processing method of the present invention. The image processing apparatus comprises an operation section 71, a 
memory 72, an input unit 73, a display unit 74, a radiotherapy treatment unit 75 and a tomographic image measurement 
apparatus 76. The operation section 71 includes a CPU, a ROM and a RAM and executes various processes such as 
a process for calculating a DRR using a color information conversion table and three-dimensional data subjected to a 
standard-deviation process. The memory 72 stores three-dimensional image data, edge-processed data, a color in- 
20 formation conversion table such as shown in Fig. 12, DRR data and the like. The input unit 73 is a device such as a 
keyboard and a mouse. The display unit 74 displays an image that is digitally reconstructed by calculation. The radi- 
otherapy treatment unit 75 is connected to the image processing apparatus through an interface section. The tomo- 
graphic image measurement apparatus 76 such as CT apparatus or NMR apparatus is connected to the image process- 
ing apparatus through an interface section. 
25 [0053] The present invention is not limited to the aforementioned embodiments. The present invention also includes 
various modifications to be made within the scope of the claims annexed to and forming a part of this specification. 



Claims 

An image processing apparatus comprising: 

a memory for storing three-dimensional data (72); 

a table (Fig. 13, 72) for allocating color information so as to correspond to voxel values ; 
an operation section (71) for generating a DRR (Digitally Reconstructed Radiograph) by interpolating voxel 
values on a ray irradiated from a radiation source and adding up color information along the ray; and 
a display unit (74) for displaying said generated DRR image; 

■ *\ 

characterized by 

said operating section (71) for adding up the color information along the ray, the color information being obtained 
by converting the interpolated voxel values based on said table. 

An image processing apparatus according to claim 1 , wherein said operation section (71 ) for generating said DRR * 
includes means for setting a plane (6) of said DRR so as to coincide with a plane including an isocenter (4). 

An image processing apparatus according to claim 2, further comprising means for expressing said three-dimen- 
sional data using a standard deviation. 

An image processing apparatus according to claim 3, wherein said operation section (71) is used for calculating 
the standard deviation of said three-dimensional data and then calculating a second version of three-dimensional 
data in which said three-dimensional data is expressed using the standard deviation. 

An image processing apparatus according to claim 2, wherein said operation section (71 ) for generating said DRR 
includes means for projecting values of voxels towards said radiation source, the voxels being located opposite 
to said radiation source with respect to the plane (6) of said DRR. 

An image processing apparatus according to claim 1 , wherein said operation section (71 ) for generating said DRR 
includes means for setting a resolution (52) of said DRR to a desired value. 



35 
. 40 

2. 

45 

3. 
4. 

50 

5. 

55 

6. 



7 



EP 0 884 697 B1 



7. An image processing apparatus according to claim 1 , wherein said operation section (71 ) for generating said DRR 
includes means for setting a sampling interval (53) between voxels on a ray connecting said radiation source to a 
pixel of said DRR to a desired value. 

5 8. An Image processing apparatus according to claim 1 , wherein said operation section (71 ) for generating said DRR 
includes 

means for partially selecting pixels of said DRR and calculating values of the selected pixels, and 
means for reconstructing said DRRs (62, 63, 64) using said partially calculated pixel values, wherein 
said display unit (74) includes means for sequentially displaying said DRRs reconstructed by using 'said 
'0 partially calculated pixel values. 

9. An image processing method/comprising: 

a step (101) of generating three-dimensional data; 
'5 a step (102) of setting radiation conditions; 

a step (103) of setting a table for allocating color information so as to correspond to voxel values ; 

a step (106) of interpolating voxel values on a ray irradiated from a radiation source; 

a step (107, 108, 104, 105) of generating a DRR (Digitally Reconstructed Radiograph) by adding up color 

information along said ray, the color information being obtained by converting the interpolated voxel values 
20 based on said table; and 

a step (109) of displaying said DRR generated. 

10. An image processing method according to claim 9, wherein said step of generating said DRR includes a step of 
setting a plane (6) of said DRR so as to coincide with a plane including an isocenter (4). 

25 

11. An image processing method according to claim 9, further comprising a step (201) of expressing said three-di- 
mensional data by calculating a standard deviation. 

12. An image processing method according to claim 11 , wherein the step of expressing said three-dimensional data 
30 comprises 

a step (201 ) of calculating a second version of three-dimensional data by calculating the standard deviation 
of said three-dimensional data. 

13. An image processing method according to claim 9, wherein said step of generating said DRR includes a step of 
35 projecting values of voxels towards said radiation source, the voxels being located opposite to said radiation source 

with respect to the plane (6) of said DRR. 

14. An Image processing method according to claim 9, wherein said step of generating said DRR includes a step of 
setting a resolution (52) of said DRR to a desired value. 

• 40 

15. An image processing method according to claim 9, wherein said step of generating said DRR includes a step 
setting a sampling interval (53) between voxels on a ray connecting said radiation source to a pixel of said DRR 
to a desired value. 

45 16. An image processing method according to claim 9, wherein said step of generating said DRR includes 
a step of partially selecting pixels of said DRR and calculating values of the selected pixels, and 
a step of reconstructing said DRRs (62, 63, 64) using said partially calculated pixel values, wherein 
said step of displaying includes a step of sequentially displaying said DRRs reconstructed by using said 
partially calculated pixel values. 

50 

Patentanspruche 

1. Bildverarbeitungsvprrichtung, mit: 

55 

einem Speicher zum Speichern dreidimensionaler Daten (72); 

airier Tabeiie (Fig. 13, 72) zum Zuordnen von Farbinformationen, so dass sie Volumenelement- bzw. Voxel- 
Werten entsprechen; 
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einer Operationssektion (71) zum Erzeugen eines DRR (digital rekonstruiertes Rontgenbildfdurch Interpolie- 
ren von Voxel-Werten auf einem von einer Strahlungsquelle eingestrahlten Strahl und Aufaddieren von Far- 
binformationen entlang dem Strahl; und 

einer Anzeigeeinheit (74) zum Anzeigen des erzeugten DRR-Bildes; 
gekennzeichnet durch 

die Operationssektion (71) zum Aufaddieren der Farbinformationen entlang dem Strahl, wobei die Farbin- 
lormationen erhalten werden, indem die interpolierten Voxel-Werte basierend auf der Tabelle umgewandelt wer- 
den. 

2. Bildverarbeitungsvorrichtung nach Anspruch 1, worin die Operationssektion (71) zum Erzeugen des DRR Einrich- 
tungen zum Festlegen einer Ebene (6) des DRR enthalt, so dass sie mit einer ein Isozentrum (4) enthaltenen 
Ebene ubereinstimmt. 

3. Bildverarbeitungsvorrichtung nach Anspruch 2, femer mit einer Einrichtung zum Ausdrucken der dreidimensiona- 
len Daten unter Verwendung einer Standardabweichung. 

4. Bildverarbeitungsvorrichtung nach Anspruch 3, worin die Operationssektion (71 ) verwendet wird, urn die Standar- 
dabweichung der dreidimensionalen Daten zu berechnen und dann eine zweite Version dreidimensionaler Daten 
zu berechnen, in der die dreidimensionalen Daten unter Verwendung der Standardabweichung ausgedruckt wer- 
den. 

5. Bildverarbeitungsvorrichtung nach Anspruch 2, worin die Operationssektion (71) zum Erzeugen des DRR eine 
Einrichtung enthalt, um Voxel-Werte in Richtung auf die Strahlungsquelle zu projizieren, wobei die Voxel der Strah- 
lungsquelle bezugiich der Ebene (6) des DRR gegenuberliegen. 

6. Bildverarbeitungsvorrichtung nach Anspruch 1, worin die Operationssektion (71) zum Erzeugen des DRR eine 
' Einrichtung zum Festlegen einer Auflosung (52) des DRR auf einen gewunschten Wert enthalt. 

7. Bildverarbeitungsvorrichtung nach Anspruch 1, worin die Operationssektion (71) zum Erzeugen des DRR eine 
Einrichtung zum Festlegen eines Abtastintervalls (53) zwischen Voxet auf einem die Strahlungsquelle mit einem 
Pixel des DRR verbindenden Strahl auf einen gewunschten Wert enthalt. 

8. Bildverarbeitungsvorrichtung nach Anspruch 1 , worin die Operationssektion (71 ) zum Erzeugen des DRR enthalt 

eine Einrichtung zum partiellen Auswahlen von Pixel des DRR und Berechnen von Werten der ausgewahlten 
Pixel und 

eine Einrichtung zum Rekonstruieren der DRRs (62, 63, 64) unter Verwendung der teilweise berechneten 
Pixelwerte, worin 

die Anzeigeeinheit (74) eine Einrichtung zum sequentiellen Anzeigen der unter Verwendung der teilweise 
berechneten Pixelwerte rekonstruierten DRRs enthalt. 

9. Bildverarbeitungsverfahren, aufweisend: 

einen Schritt (101) zum Erzeugen dreidimensionaler Daten; 
einen Schritt (102) zum Festlegen von Strahlungsbedingungen; 

einen Schritt (103) zum Festlegen einer Tabelle zum Zuordnen von Farbinformationen, so dass sie Voxel- 
Werten entsprechen; 

einen Schritt (106) zum Interpolieren von Voxel-Werten auf einem von einer Strahlungsquelle eingestrahlten 
Strahl; 

einen Schritt (107, 108, 104, 105) zum Erzeugen eines DRR (digital konstruierten Rontgenbild), indem Farb- 
informationen entlang dem Strahl aufaddiert werden, wobei die Farbinformationen erhalten werden, indem 
die interpolierten Voxet-Werte basierend auf der Tabelle umgewandelt werden; und 
einen Schritt (109) zum Anzeigen des erzeugten DRR. 

10. Bildverarbeitungsverfahren nach Anspruch 9, worin der Schritt zum Erzeugen des DRR einen Schritt zum Festle- 
gen einer Ebene (6) des DRR enthalt, so dass sie mit einer ein Isozentrum (4) enthaltenden Ebene ubereinstimmt. 

1 1 . Biidverarbeitungsverfahren nach Anspruch 9, ferner mit einem Schritt (201 ) zum Ausdrucken der dreidimensiona- 
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len Daten durch Berechnen einer Standardabweichung. 

12. Bildverarbeitungsverfahren nach Anspruch 11, worin der Schritt zum Ausdrucken der dreidimensionalen Daten 
aufweist 

einen Schritt (201) zum Berechnen einer zweiten Version dreidimensionaler Daten durch Berechnen der 
Standardabweichung der dreidimensionalen Daten. 

13. Bildverarbeitungsverfahren nach Anspruch 9, worin der Schritt zum Erzeugen des DRR einen Schritt zum Proji- 
zieren von Voxel-Werten in Richtung aut die Strahlungsquelle enthalt, wobei die Voxel der Strahlungsquelle be- 
zuglich der Ebene (6) des DRR gegeniiberliegen. 

14. Bildverarbeitungsverfahren nach Anspruch 9, worin der Schritt zum Erzeugen des DRR einen Schritt zum Festle- . 
gen einer Auflosung (52) des DRR auf einen gewiinschten Wert enthalt. 

15. Bildverarbeitungsverfahren nach Anspruch 9, worin der Schritt zum Erzeugen des DRR einen Schritt einschlieBt, 
der ein Abtastintervall (53) zwischen Voxel auf einem die Strahlungsquelle mit einem Pixel des DRR verbindenden 
Strahl auf einen gewunschten Wert festlegt. 

16. Bildverarbeitungsverfahren nach Anspruch 9, worin der Schritt zum Erzeugen des DRR einschlieBt 

einen Schritt zum teilweisen Auswahlen von Pixel des DRR und Berechnen von Werten der ausgewahlten 
Pixel und 

einen Schritt zum Rekonstruieren der DRRs (62, 63, 64) unter Verwendung der teilweise berechneten Pi- 
xelwerte, worin 

der Schritt zum Anzeigen einen Schritt zum sequentiellen Anzeigen der unter Verwendung der teilweise 
berechneten Pixeiwerte rekonstruierten DRRs einschlieBt. 



Revendications 

1. Dispositif de traitement d'images comportant : 

une m6moire pour m£moriser de donnees tridimensionnelles (72), 

une table (figure 13, 72) pour affecter des informations de couleur afin qu'elles correspondent a des valeurs 
de voxels, 

une section operationnelle (71) pour generer une Radiographie NumSriquement Reconstitute (DRR) en in- 
terpolant des valeurs de voxels sur un rayon irradie par une source de rayonnement et en additionnant des 
informations de couleur le long du rayon, et 
une unite d'affichage (74) pour afficher ladite image DRR g6n£r£e, 

caracterise par 

ladite section operationnelle (71) destinee a addittonner les informations de couleur le long du rayon, les 
informations de couleur etant obtenues en convertissant les valeurs de voxels interpolees sur la base de ladite 
table. 

2. Dispositif de traitement d'images selon la revendication 1, dans lequel la section operationnelle (71) destinee a 
g 6n£rer ladite DRR inclut des moyens pour Stablir un plan (6) de ladite DRR afin qu'il coincide avec un plan incluant 
un isocentre (4). 

3. Dispositif de traitement d'images selon la revendication 2, comportant de plus des moyens pour exprimer lesdites 
donnees tridimensionnelles en utilisant un ecart type. 

4. Dispositif de traitement d'images selon la revendication 3, dans lequel ladite section operationnelle (71 ) est utilisee 
pour calculer I'ecart type desdites donnees tridimensionnelles et calculer ensuite une seconde version de donnees 
tridimensionnelles dans laquelle lesdites donn6es tridimensionnelles sont exprimees en utilisant I'ecart type, 

5. Dispositif de traitement d'images seion la revendication 2, dans lequel ladite section operationnelle (71 ) destinee 
a generer iadite DRR inciut des moyens pour projeier des vaieurs de voxeis vers iadiie source de rayonnemeni, 
les voxels etant situ6s a ('oppose de ladite source de rayonnement par rapport au plan (6) de ladite DRR. 
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, 6. Dispositif de traitement d'images seton la revendication 1 , dans lequel ladite section operationnelle (71 ) destin6e 
a gen6rer ladite DRR inclut des moyens pour etablir une resolution (52) de ladite DRR a une valeur voulue. 

7. Dispositif de traitement d'images selon la revendication 1 , dans lequel ladite section operationnelle (71 ) destin6e 
5 a g6n6rer ladite DRR inclut des moyens pour etablir a une valeur voulue un intervalle d'echantillonnage (53) entre 

des voxels sur un rayon reliant ladite source de rayonnement a un pixel de ladite DRR. 

8. Dispositif de traitement d'images selon la revendication 1, dans lequel ladite section operationnelle (71) destinee 
a generer ladite DRR inclut 

w des moyens pour seiectionner partiellement des pixels de ladite DRR et calculer des valeurs des pixels 

s6lecttonnes, et 

des moyens pour reconstituer lesdites DRR (62, 63 ( 64) en utilisant lesdites valeurs de pixels partiellement 
calculees, dans tequel 

ladite unite d'afiichage (74) inclut des moyens pour afficher sequentiellement lesdites DRR reconstitutes en 
is utilisant lesdites valeurs de pixels partiellement calculees. 

9. Procede de traitement d'images, comportant : 

une etape (101) de generation de donnees tridimensionnelles, 
20 une etape (1 02) d'etablissement de conditions de rayonnement, 

une etape ( 1 03) d'etablissement d'une table pour affecter des informations de couleur af in qu'elles correspon- 
dent a des valeurs de voxels, 

une etape (106) ^interpolation des valeurs de voxels sur un rayon irradie par une source de rayonnement, 
une etape (107, 108, 104, 105) de generation d'une DRR (Radiographie Numeriquement Reconstitute) en 
25 additionnant des informations de couleur le long dudit rayon, les informations de couleur etant obtenues en 

convertissant les valeurs de voxels interposes sur la base de ladite table, et 
une etape (109) d'affichage de ladite DRR gentrte. 

10. Procede de traitement d'images selon la revendication 9. dans lequel ladite etape de generation de ladite DRR 
30 inclut une etape d'etablissement d'un plan (6) de ladite DRR af in qu'il coincide avec un plan incluant un isocentre (4). 

11. Procede de traitement d'images selon la revendication 9, comportant de plus une etape (201 ) consistant k exprimer 
lesdites donntes tridimensionnelles en calculant un ecart type, 

35 12. Procede de traitement d'images selon la revendication 11 , dans lequel retape consistant a exprimer lesdites don- 
nees tridimensionnelles comporte 

une etape (201) de calcul d'une seconde version de donnees tridimensionnelles en calculant I'ecart type 
desdites donnees tridimensionnelles. 

to 13. Procede de traitement d'images seton la revendication 9, dans lequel ladite etape de generation de ladite DRR 
inclut une etape de projection de valeurs de voxels vers ladite source de rayonnement, les voxels etant situes k 
I'oppose de ladite source de rayonnement par rapport au plan (6) de ladite DRR. 

14. Procede de traitement d'images selon la revendication 9, dans lequel ladite etape de generation de ladite DRR 
45 inclut une etape d'etablissement d'une resolution (52) de ladite DRR a une valeur voulue. 

15. Procede de traitement d'images selon la revendication 9, dans tequel ladite etape de generation de ladite DRR 
inclut une etape d'etablissement a une valeur voulue d'un intervalle d'echantillonnage (53) entre des voxels sur 
un rayon reliant ladite source de rayonnement a un pixel de ladite DRR. 

50 

16. Procede de traitement d'images selon la revendication 9, dans lequel ladite etape de generation de ladite DRR 
inclut 

une etape de selection partielle des pixels de ladite DRR et de calcul des valeurs desdits pixels selectionnes, 

et 

55 une etape de reconstitution desdites DRR (62, 63, 64) en utilisant lesdites valeurs de pixels partiellement 

calculees, dans tequel 

ladite etape d'auichaye inciui une eiape d : affichage sequenttei desdites DRR reconstitutes en utilisant les- 
dites valeurs de pixels partiellement calculees. 
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FIG. 13 
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